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Abstract 

We consider the scattering by a one-dimensional random potential and derive the 
probability distribution of the corresponding Wigner time delay. It is shown that the 
limiting distribution is the same for two different models and coincides with the one 
predicted by random matrix theory. It is also shown that the corresponding stochastic 
process is given by an exponential functional of the potential. 
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1 Introduction 



The concept of time delay introduced long ago by Eisenbud and Wigner has recently 
received renewed attention. In a scattering process the time delay r(k) is related to the 
time spent in the interaction region by a wavepacketQ with energy peaked at E = k 2 . It can 
be expressed in terms of the derivative of the S matrix with respect to the energy. If the 
system is described by a finite number N of channels, the N time delays are the eigenvalues 
of the matrix — iS^j^. In the context of chaotic scattering a statistical approach based on 
random matrix theory allows a determination of the complete distribution of time delays 
(for a review see Q). This problem was first studied in [|j by a super symmetric approach, 
and in || by using a statistical analysis. This latter work provided a physical derivation 
of the one-channel case for (3 = 2 and gave the distribution for other universality classes. 
It then recently served as a starting point for || where the iV-channel distribution is 
shown to be given by the Laguerre ensemble of random matrix theory. A slightly different 
approach has been developed in the context of quasi one-dimensional mesoscopic systems. 
Quite recently a resonant transmission model was proposed to account for the different 
behaviours observed in the metallic and insulating regimes |5J. 

In this paper we consider a one-dimensional model with a random potential whose 
support is a finite segment of length L. A hard wall condition at the origin reduces the 
problem to a scattering problem on the half line. Using standard methods we write down a 
set of two coupled stochastic differential equations satisfied by the phase and its derivative 
with respect to the energy which is related to the time delay. In the high energy (or weak 
disorder) limit we may take advantage of the fact that the phase is a rapid variable with 
a uniform distribution. One finally ends up with a stochastic differential equation for the 
time delay which coincides with the one derived in [^, ||] and which can be integrated 
explicitly in the case where the potential is white noise. The resulting expression is given 
in terms of an exponential functional of a Brownian motion. An analysis of the same 
model with a different type of disorder yields the same functional and consequently the 
same stationary distribution. Our work therefore brings some support to the conjecture 
that these distributions are universal. Another aim of this work is to draw attention to 
the fact that the same probability distribution also arises in other physical problems, for 
instance in the context of diffusion in a random medium. 



2 The model with a white noise potential 

Consider the Schrodinger equation on the half line x > 

d 2 

— 7ju(x) + V(x)u(x) = k 2 u(x) . (1) 

dx z 

We assume that the potential has its support on the interval [0, L] and impose Dirichlet 
boundary condition at the origin u(0) = 0. Therefore for x > L scattering states of the 
form u(x) = e~ zkx + e ^x+ie(k) re p resen ^ th e superposition of an incoming plane wave 
and a reflected plane wave. In this case the reflection coefficient r(k) = e ld ^) is of unit 
modulus since there is only backward scattering. All information on the scattering process 

4 For a detailed and critical analysis of this concept see jl]] and [[|. 
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is contained in the phase shift 9{k). In particular the Wigner time delay is given by 

Our aim is to find the probability distribution of r(k) in the case where the potential is a 
Gaussian white noise such that 

(V(x)) = 

(V(x)V(y)) = aj(x-y). (3) 

A widely used method to deal with such problems is the invariant embedding approach 
flO| , ^ . Here we will use a more straightforward derivation following [11]. 

A 

dx 

equation 

V(x) - k 2 -z{xf . (4) 



The Ricatti variable z(x) = ■4-\n\u(x)\ satisfies the first order nonlinear differential 

dz ( X ) s~, , ,2 ,2 



dx 

In the region where V{x) vanishes one has 

*\ x > ~ i i e 2ikx+ie(k) ■ 

Inverting this relation we may express the reflection coefficient in terms of the Ricatti 
variable 

2ikx+i9(k) _ ^ + z ( x ) (r \ 

~ ik- z(x) • [b) 

Although this relation only holds in the region x > L, we may take it as a definition of 
the variable i/j(x) 

Q iip(x) _ + z ( x ) fj\ 
ik — z(x) 

in the region < x < L. By using equation (Q) one gets the differential equation satisfied 

by i>(x) 

^tM. = 2k- r(l + cos Mx))V(x) . (8) 
dx k V / 

Since the Ricatti variable is continuous even if the potential has a discontinuity this implies 
that the phase shift 6{k) is given by if)(L, k) — 2kL. The reflection condition at the origin 
gives the boundary condition -0(0) = vr. 

Let us introduce the derivative of ip with respect to k: Z(x, k) = -^^{x, k) in terms of 
which the time delay reads 

rW = ^|^ (9) 

The differential equation for Z now reads 



dZ(x) 
dx 



— « ( 1 + cos ip(x)) + T Z{x) sin ip(x) 
k z V / k 



V(x) . (10) 



From the boundary condition at x = we may set Z(0) = 0. 



3 



If V(x) is a Gaussian white noise equations (|8 , 10 ) are coupled stochastic differential 
equations (in the Stratonovich sense) which give the phase and the time delay. 

One then may pass from these two stochastic differential equations to a Fokker-Planck 
equation for the probability density of ip and Z. In the high energy limit k 3> cr*/ 3 one can 
show that the variable ip is a rapid variable uniformly distributed on the interval [0, 2tt]. 
Moreover, since the rapid variable ip and the slow variable Z decorrelate in this limit, then 
one may average over the rapid variable and eventually get the following Fokker-Planck 
equation for Z (a more detailed description of the procedure will be given in section ||) 

dP(Z;x) _ d 
~ ~dZ 



dx 



2)P(Z:x) 



+ WWz z Wz ZP ^ 



(ii) 



Up to an inessential termQ, this equation coincides with the one derived in § and g. It 
belongs to a more general class of Fokker-Planck equations that appear in the study of 
exponential functionals of Brownian motion with drift. 



A representation of r as an exponential functional of the 
Brownian motion 



Equation (11) describes a special case of a class of stochastic processes which have been 
studied extensively in the mathematical [12| as well as in the physical literature. It can be 
cast into the equation studied by Schenzle et al. in the context of multiplicative stochastic 
processes pM. M ore recently it was shown to arise in the context of diffusion in a random 
medium |14|, [Dj]. The distribution of the flux of particles [16, |T^] in a disordered sample of 
finite length or the waiting time distribution may be obtained by solving a generalization 
of equation (11). From our previous work [17, 18] one may write the general solution of 
(|ll| ) in the form 



7-./ T\ ^ A 2 _A_ f°° , ifliJl S . TTS 

Z z ttZ Jo 1 + s z 2 



W, is - 



(12) 



where W^ v (z) is a Whittaker's function and A = is the localization length at high 

energy (k S> cr}^ 3 )- In the limit L — > cxd the first term gives the limiting distribution 
which, as announced, coincides with the one channel distribution obtained in 0. It is also 
in agreement with the result of Jayannavar et al. @ in the high energy limit. 

The general solution given above allows to study finite size effects. The first correction 
to the stationary distribution for Z may be calculated from (|i~2|) 



P(Z-L) ~ A e - 

L>A Z z 



A /7T /A\ 3 / 2 



2 \L 



A. /A 
2A Wxjq ( - 



A 

e 2z 
Z 



(13) 



One can also compute all the moments of the variable Z (or r) |17], 18 1. For n greater 



than one they all diverge exponentially as a function of the length of the sample. 

(Z(L)) = 2L (14) 



{Z(LY 



(n - 2)1 
(2n - 2)! 



A n e 2n(n-1)A 



(15) 



we do not write this term because it is of the same order as terms we neglected in the approximation 
of the decorrelation of ip and Z. 
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It is interesting to remark that the limiting distribution is not of the log-normal form as 
could have been expected |J from a simple resummation of the most divergent part of the 
moments. 

Another by-product of our earlier work is to provide a representation of the process 
to which the random variable r obeys. The corresponding stochastic differential equation 
associated with ( |Tl|) may be written in the following form 



d ^= 2 -w^ + irk m:) (16) 

where V(x) is the white noise (|3|). By integration one obtains Z and from (||) the time 
delay 




T (k) = - dx\e Jx ik 'J -1 . (17) 



A more satisfactory procedure to derive this formula is to start from the stochastic dif- 
ferential equations (^[H]) and construct r without using a Fokker-Planck equation. This 
approach was used by Faris and Tsay |l9| . 

When writing this expression for r(k) it should be kept in mind that although (|l~7| ) is 
not true for every realization of V(x), it nevertheless captures all the statistical properties 
of the process. In particular it allows the study of the various limiting cases. Since the 
potential V(x) is a white noise, the first term in the exponential, which is a Brownian 



motion, is typically of order ya g L/k 2 . The argument of the exponential is roughly a 
function of the ratio of the localization length X(k) = 8k 2 /a g and the size of the system L. 
In particular if A <C L the time delay is approximatively equal to —L/k which is twice the 
length of the sample divided by the speed of the particle. This value for the time delay 
corresponds to a reflection at the sample edge. This agrees with the fact that if A <C L, 
one expects that the particle will only enter the sample for a negligeable length compared 
to L. In the other regime A > L expression (17) shows that the time delay is roughly 



zero. This is consistent with the fact that in this regime the sample is almost transparent 
to the particle. 

Moreover if one considers the regime where A <C L, the stationary distribution for 
Z(V) gives a typical value Z tjp ~ 2A corresponding to a typical time delay r typ ~ — -^p. 
Such a value for the time delay means that the particle enters in the disordered region for 
a typical length A. 

One may wonder to what extent the representation of Z as an exponential functional of 
the potential depends on the specific form of the noise that we have considered. Although 
we are not able to give a definite answer to this question we nevertheless notice that for 
any realization of V(x) one may integrate equation ( |io|) and get 

Z(L) = 2 £ dx (l+ ] ^{l + cos^x))^e^^ dx ' V{x ' )s ^ x ' ) (18) 
(A;) * \ ("dx Uj>' v W™^-l\ . (19) 



hence 

T 
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Up to a constant drift, at high energy, this is essentially of the form given in (|17|). The 
distribution and the asymptotic behaviour of exponential functionals of more general pro- 
cesses like Levy processes has been obtained in the mathematical literature [20, 21 1. There 
it is shown that a Poisson process gives a limiting distribution whose tail still decays al- 
gebraically. 



4 The supersymmetric model 
4.1 The model 

In this section we consider a different model for which one can also compute exactly the 
time delay distribution in the weak disorder limit. We show below that r(fc) is distributed 
with the same law as in the previous case. 

We consider the one-dimensional Schrodinger Hamiltonian 

H = -^ + ^{x)+^{x). (20) 
ax z 



As explained at length in |T4], g^, 23, [24J] this model arises in diverse areas of quantum 
mechanics ranging from the study of solitons in polymers to the study of classical diffusion 
in a random medium. Recently it was used in the context of one-dimensional spin chains 
[25] and isospectral periodic potentials [26]. 

In the following we consider the case where (j)(x) is a Gaussian white noise with the 
moments 

<<K*)> = 

(<Kx)<l>(v)) = aS(x-y). (21) 
The supersymmetric Hamiltonian (p0|) can be rewritten in the factorized form 

H = Q ] Q (22) 

where 

Q = -£ + c(>(x) (23) 
ax 

Qt = ± + 4>(x). (24) 
ax 

The density of states and the localization length for this model were first obtained in 



[27] and then rediscovered independently in [28]. In contrast with the previous model, in 
the high energy limit (k 3> a) the localization length reaches the constant value X(k) = -. 

The two other length scales of the problem are the size of the disorder region L and 
the de Broglie wavelength k^ 1 . In the following we choose to work in a high energy limit, 
i.e. when k ^> a and k ^> L _1 . 

As in the previous section, the potential is non zero on the interval [0, L\ and there is 
a reflection condition at the origin. 
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4.2 Stochastic differential equations for the phase and the time delay 

Using the factorization of the supersymmetric Hamiltonian, one may decouple the sta- 
tionary Schrodinger equation {^—-jj^z + 4> 2 {x) + (j>'(x)j u(x) = k 2 u(x) into two first order 
equations: 

dU ^ ■= 4>(x)u{x) - kv{x) (25) 



dx 
dv(x) 

dx 



~4>(x)v(x) + ku{x) . (26) 



If <p{x) has a discontinuity the two functions u(x) and v(x) are continuous. This suggests 
the introduction of the Ricatti variable ((x) = It is straightforward to see that C( x ) 
obeys the following first order non linear differential equation: 

= k _ 2<P(x)C(x) + k ((x) 2 . (27) 

In the region where (j)(x) is vanishing (x > L), the stationary scattering states can be 
expressed as 

u(x) = e~ ikx + e ikx+i6 ^ (28) 

where 9{k) is the phase shift. 

This leads to the change of variable: 

i p,2ikx+ia{x) 

C(x) = i 17- l s (29) 

with the phase shift 9{k) given by a(L, k). Instead of a{x) it is in fact more convenient 
to introduce the variable ift(x) = a(x) + 2kx in terms of which the Ricatti variable is 
C(x) = tan(^(x)/2). Equation (p?|) then gives 

= 2k-2(j)(x)smi>{x) . (30) 

dx 

This gives the evolution of the phase and consequently the distribution of the nodes of the 
wave function u(x) from which one can get the density of states (compare with equation 
(2.5) of H). 

As in the previous section one may introduce Z(x,k), the derivative of if)(x,k) with 
respect to k, which satisfies 

= 2 - 2^(x) Z(x) cos if;(x) (31) 

dx 

again the initial conditions for the two variables are ip(0) = it and Z(0) = 0. 

Integrating the coupled equations (3C,31) between and L gives the phase shift and 
the time delay. 
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4.3 The phase distribution 

At this stage the formalism developed to compute the phase shift and the time delay is 
quite general. We now consider the case where (f>(x) is white noise. Equation (^) is then 
a stochastic differential equation in the Stratonovich sense from which we can write down 
a Fokker-Planck equation for the probability density of ip : 

d - E ^ 1 = -2k^P(^x)+2a^sm^sm^-,x) . (32) 
The stationary distribution for ip is given by: 

PsW=M ■ , / (33) 
sin tp Jip sm yy 

where N is a normalization constant related to the integrated density of states per unit 
length N{E) = 2kP s (ir) = 2oM (see equation (2.18) of |22| and reference therein). One 
can extract from (^) the high energy expansion of the stationary distribution: 

PsW=A/ "i{ 1+ S sin2 ^ + SG sin2 2lp - sin4 ^) + (i) } ■ (34) 

This shows that in the high energy limit the phase is a uniformly distributed variable on 
[0, 2tt]. This could have been guessed directly from the differential equation for tp since in 
equation ( |30|) one expects that for k ^> a the first term will dominate. 

4.4 The time delay 

Since if) is a rapidly varying variable, one expects that the two coupled differential equations 
will reduce to only one, after averaging over the rapid variable. For this purpose, the next 
natural assumption is the decorrelation of the two variables i/j and Z in the high energy 
limit. We first integrate perturbatively equations (^ 31) to show that this is indeed the 
case: 

i/,(x) = vr + 2kx + 2 / dx' (j){x') sin 2kx + • • • (35) 



Since the integral is of the order ^/ax, this expansion is valid if \fox <C 1 . The computation 
of the autocorrelation function { r ip{x) r ip{y)) — {il){x)){il){y)) then shows that the variable tp 
behaves in this limit as a Brownian motion with a drift: 

tp[x) ~ vr + 2kx + y/2 f dx' (p(x') . (36) 
Jo 

One can perform the same approximation for Z and compute the correlation function with 
the help of the two expansions 

(ip(x)Z(x)) — {ip(x)){Z(x)) ~ — sin Akx — ax cos Akx^j . (37) 

From which we deduce that the correlations between if) and Z will vanish in the high 
energy limit. 
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From the two stochastic differential equations ( P0|j31| ) we may write the Fokker-Planck 
equation for the joint probability density P(tp, Z; x): 

dP dP dP [ d d d d 

—— = —2k— 2— \-2a { — — svaip—— svaip P + — — sin^— — Zcos ip P 

ox dip oZ ( dip dip dip dZ 

d d d d 1 

+—Z cos ip— sinip P + —Z cos ip—Z cos ip Pj . (38) 

Using the fact that ip is uniformly distributed and decorrelated from Z in the high energy 
limit, we may average over ip since the joint probability density factorises as P(ip, Z; x) ~ 
■£-P(Z;x). This leads to the following Fokker-Planck equation for P(Z;x): 

^> = §z ^ z -«"» + ^4 XPR,) (39) 



which is the same as equation (11) provided one replaces by ^. In particular this 
implies that the distribution for the variable Z is still given by equation ( |l2|) provided A 
is the localization length for the super symmetric model. 

One may thus give a representation for the time delay, valid for both models: 

7 kj ( e 1 ) ( °) 

where r](x) is the white noise of variance 1 entering in the potential (equal to -^-V{x) in 
the first case and -^cp(x) in the second) and A is the localization length. 

5 Conclusion 

We have shown in this paper that, for two models of random potential, the time delay 
exhibits the same distribution. These statistical properties are in agreement with other 
approaches, underlying the universality of such properties. 

We have also demonstrated that in both cases the time delay may be expressed as an 
exponential functional of the white noise which enters in the potential. 

It would be interesting to extend this approach to the multichannel case for which 
the stationary distribution has been recently obtained Q. Another point that deserves 
attention is the fact that the resulting distribution for the time delay is exactly of the 
same form as the waiting time distribution that occurs in the context of classical diffusion 



in a random potential [14|. The supersymmetric model may help to explore this relation. 
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